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The vulcanization transition — the crosslink-density-controlled equilibrium phase transition from 
the liquid to the amorphous solid state — is explored analytically from a renormalization group 
perspective. The analysis centers on a minimal model which has previously been shown to yield 
a rich and informative picture of vulcanized matter at the mean-field level, including a connection 
with mean-field percolation theory (i.e. random graph theory). This minimal model accounts for 
both the thermal motion of the constituents and the quenched random constraints imposed on their 
motion by the crosslinks, as well as particle-particle repulsion which suppresses density fluctuations 
and plays a pivotal role in determining the symmetry structure (and hence properties) of the model. 
A correlation function involving fluctuations of the amorphous solid order parameter, the behavior 
of which signals the vulcanization transition, is examined, its physical meaning is elucidated, and 
the associated susceptibility is constructed and analyzed. A Ginzburg criterion for the width (in 
crosslink density) of the critical region is derived and is found to be consistent with a prediction 
due to de Gennes. Inter alia, this criterion indicates that the upper critical dimension for the 
vulcanization transition is six. Certain universal critical exponents characterizing the vulcanization 
transition are computed, to lowest nontrivial order, within the framework of an expansion around 
the upper critical dimension. This expansion shows that the connection between vulcanization and 
percolation extends beyond mean-field theory, surviving the incorporation of fiuctuations in the sense 
that pairs of physically analogous quantities (one percolation-related and one vulcanization-related) 
are found to be governed by identical critical exponents, at least to first order in the departure 
from the upper critical dimension (and presumably beyond). The relationship between the present 
approach to vulcanized matter and other approaches, such as those based on gelation/percolation 
ideas, is explored in the light of this connection. To conclude, some expectations for how the 
vulcanization transition is realized in two dimensions, developed with H. E. Castillo, are discussed. 
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I. INTRODUCTION 



Whilst a rather detailed description of the vulcanization transition has emerged over the past few years within 
the context of a mean-field approximation 0-^], the picture of this transition beyond the mean- field level is less 
certain. The purpose of the present Paper is to provide a description of the vulcanization transition beyond the 
mean-field approximation via the application of renormalization group (RG) ideas to a model that incorporates both 
the quenched randomness (central to systems undergoing the vulcanization transition) and the thermal fluctuations of 
the constituents (whose change in character is the fundamental hallmark of the transition). Our aim is to shed some 
light on certain universal properties of the vulcanization transition within the framework of the well-controlled and 
systematically improvable approximation scheme that the RG provides, viz., an expansion about an upper critical 
dimension that we shall see takes the value six. 

We remind the reader that the vulcanization transition is an equilibrium phase transition from a liquid state of 
matter to an amorphous solid state. (In addition to the technical reports cited above we refer the reader to some 

informal accounts of the physics of the vulcanization transition |^J^,^.) The transition occurs when a sufficient density 
of permanent random constraints (e.g. chemical crosslinks) — the quenched randomness — are introduced to connect the 
constituents (e.g. macromolecules) , whose locations are the thermally fluctuating variables. In the resulting amorphous 
solid state, the thermal motion of (at least a fraction of) the constituents of the liquid undergo a qualitative change: 
no longer wandering throughout the container, they are instead localized in space at random positions about which 
they execute thermal (i.e. Brownian) motion characterized by random r.m.s. displacements. 

Our approach to the vulcanization transition is based on a minimal Landau- Wilson effective Hamiltonian that 
describes the energetics of various order-parameter-field configurations, the order parameter in question having been 
crafted to detect and diagnose amorphous solidification. This order parameter and effective Hamiltonian can be derived 
(along with specific values for the coefficients of the terms in the effective Hamiltonian) via the application of replica 
statistical mechanics to a specific semi-microscopic model of randomly crosslinked macromolecular systems (RCMSs) , 
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viz., the Deam-Edwards model this procedure is described in detail in Ref. More generally, the form of the 
minimal model can be determined from the nature of the order parameter, especially its transformation properties 
and certain symmetries that the effective Hamiltonian need possess, along with the assumptions of the analyticity 
of the effective Hamiltonian and the continuity of the transition. This system-nonspecific strategy for determining 
the minimal model was applied in Rcf. JlOt . There, it was shown that by regarding the effective Hamiltonian as a 
Landau free energy one could recover from it the mean-field description of both the liquid and emergent amorphous 
solid states known earlier from the analysis of various semi- microscopic models [^,|[^ 12 . The mean- field value 



of order parameter in the solid state encodes a function rather than a number, and it possesses a certain mean- 
field "universality" by which we mean that (as the transition is approached from the amorphous solid side) both 
the exponent governing the vanishing of the fraction of constituents localized (i.e. the gel fraction) and the scaled 
distribution of localization lengths of the localized constituents turn out to depend not on the coefhcients in the Landau 
free energy but only on its qualitative structure. Support for this mean-field picture of the amorphous solid state, in 
the form of results for the localized fraction and scaled distribution of localization lengths, has emerged from extensive 
molecular dynamics computer simulations of three-dimensional, off-lattice, interacting, macromolecular systems, due 
to Barsky and Plischke |l3|Jl4|| . In order to provide a unified theory of the vulcanization transition that encompasses 
the liquid, critical and random solid states, we shall in the present work be adopting this Landau free energy as the 
appropriate Landau- Wilson effective Hamiltonian. 

We shall be focusing on the liquid and critical states, rather than the amorphous solid state, and shall therefore be 
concerned with the order-parameter correlator rather than its mean value. Along the way, we shall therefore discuss 
the physical content of this correlator, why it signals the approaching amorphous solid state, and how it gives rise to 
an associated susceptibility whose divergence mark the vulcanization transition. 

Given the apparent precision of the picture of the amorphous solid state resulting from the mean-field approxima- 
tion §,|jlO[|l|,|l|, the reader may question the wisdom of our embarking on program that seeks to go beyond the 
mean-field approximation by incorporating the effects of fluctuations. We therefore now pause to explain what has 
motivated this program. 

(i) Below six spatial dimensions, mean-field theory necessarily breaks down sufficiently close to the vulcanization tran- 
sition. Although, as we shall also see, the region of crosslink densities within which fluctuations play an important 
role is narrower for dimensions closer to (but below) six and for longer macromolecules, it is by no means necessary 
for this region to be narrow for shorter macromolecules and for lower-dimensional systems; thus, systems for which 
the fluctuation-dominated regime is observably wide certainly exist. 

(ii) Whilst there have been many successful treatments of critical phenomena beyond the mean-field approximation in 
systems with quenched randomness, these have, by and large, been for systems in which the emergent order was not 
of the essentially random type under consideration here or in the spin glass setting |l^] . Instead the emergent order 
has typically been of the type arising in pure systems, albeit perturbed by the quenched disorder. We are motivated 
here by the challenge of going beyond mean-fleld theory in the context of a transition to a structurally random state 
of matter. 

(iii) The vulcanization transition has often been addressed from the perspective of gelation/percolation theo- 
ries |l^-|20|]. Whilst this perspective can be (and certainly has been) taken beyond the mean-fleld level, it possesses 
but a single ensemble, and therefore does not incorporate the effects of both quenched randomness and thermal fluc- 
tuations ]2^ . Given that an essential aspect of the vulcanization transition is the impact of the quenched random 
constraints on the thermal motion of the constituents, the a priori identification of the vulcanization transition with 
gelation/percolation is thus a nontrivial matter. By contrast with the gelation/percolation- type of approaches, the 
analysis given in the present Paper applies directly to the vulcanization transition exhibited by thermally fluctuating 
systems and driven by quenched random constraints. It should therefore shed some light on the relevance of the 
gelation/percolation- type perspective for the vulcanization transition, as we shall discuss in Sec. VI. 

This Paper is organized as follows. In Sec. || we give a brief account of the order parameter for the vulcanization 
transition, and of the Deam-Edwards replica approach to vulcanized matter and its field-theoretic representation. 



together with a minimal field-theoretic model for the vulcanization transition. In Sec. [II we summarize the mean- 



field-level picture of the vulcanization transition, along with the picture of the amorphous solid state that emerges 



from it. In Sec. [V we discuss the order-parameter correlator and susceptibility for the vulcanization transition, 
and examine their physical content. Having established this preparatory framework, we embark, in Sec. on the 
analysis of the vulcanization transition beyond mean-field theory. We begin by examining the self-consistency of mean- 
field theory by estimating the impact of fluctuations perturbatively, which results in the construction of a Ginzburg 
criterion and the identification of six as being the appropriate upper critical dimension. We then apply a momentum- 
shell KG scheme to the minimal model, thus obtaining certain universal critical exponents in an expansion around 
six dimensions. Finally, in Sec. VI we give some concluding remarks in which we discuss connections between our 
approach and those based on gelation/percolation, and we examine the role played by thermal fluctuations, especially 
in lower spatial dimensionalities. In three appendices we provide technical details associated with the derivation of 



2 



the Ginzburg criterion, we investigate the effects of various fields and vertices omitted from the minimal model, and 
we present the full derivation of the RG flow equations. 



II. MODELING THE VULCANIZATION TRANSITION 



The purpose of the present section is to collect together the basic ingredients of our approach to the vulcanization 
transition, including the order parameter, underlying semi-microscopic model, replica field theory, and minimal model. 
All these elements have been discussed in detail elsewhere, and we shall therefore be brief. As the reader will see, 
although its construction follows a quite conventional path, the theory does possess some intricacies. We shall therefore 
take various opportunities to shed some light on the physical meaning of its various ingredients. 

Although most of our results are not specific to any particular system undergoing a vulcanization transition, in 
order to make our presentation concrete we shall discuss the physical content for, and use notation specific to, the case 
of RCMSs. We shall follow closely the notation of Ref. [Q and, accordingly, we shall adopt units of length in whi ch 
the characteristic size of the macromolecules is unity (except in our discussion of the Ginzburg criterion, Sec. VA). 



A. Order parameter for the vulcanization transition 



The appropriate order parameter for the vulcanization transition, capable inter alia of distinguishing between the 
liquid and amorphous solid states, is the following function of A wavevectors {k^, k^, • • • , k'^}: 



[n^^ dsiexpik^ ■Cj{s))^{expik^ ■Cj{s))^---{expik' 



(2.1) 



where N is the total number of macromolecules, Cj(s) (with j = 1,...,7V and < s < 1) is the position in d- 
dimensional space of the monomer at fractional arclength s along the j"^ macromolecule, {■ ■ ■)x denotes a thermal 
average for a particular realization x of the quenched disorder (i.e. the crosslinking), and [• • •] represents a suitable 
averaging over this quenched disorder. It is worth emphasizing that the disorder resides in the specification of what 
monomers are crosslinked together: the resulting constraints do not explicitly break the translational symmetry of 
the system. In the liquid state, for each monomer {j, s) the thermal average (expzk • Cj(s))^ takes the value S^^q and 

thus the order parameter is simply Ila=i '^k° o- other hand, in the amorphous solid state we expect a nonzero 

fraction of the monomers to be localized, and for such monomers (exp i'k-Cj{s))x takes the form ,5)(k) exp ik-bj (s), 
i.e., a random phase- factor determined by the random mean position hj{s) of the monomer (j, s) times a random 
Debye- Waller factor p(j^s) (k) describing the random extent to which the monomer is localized. As reviewed in Sec. 3 of 
Rcf . 101 , by choosing the wavevectors {k" }f^to satisfy the constraint k^ + k^ + • • • + k"^ = the random phase- factors 

and hence the order parameter is capable of distinguishing between 



are eliminated from the order parameter (2.1 



the liquid and amorphous solid states and, furthermore, characterizing the randomness of the localization through its 
dependence on the collection of wavevectors. 



B. Replicated semi-microscopic model of vulcanized macromolecular systems 

Following Deam and Edwards by (i) starting from a semi- microscopic Hamiltonian describing a system of 
macromolecules interacting via an excluded-volume interaction, (ii) introducing the random constraints imposed by 
crosslinking, and (iii) averaging over the quenched disorder using the replica technique (with a physical choice for 
the distribution of the disorder which leads to an additional replica), one arrives at the disorder-averaged, replicated 
partition function (for details, see Sec. 4 of Ref. Q) 

[Z"](x(exp(-HP+i))^^^ , (2.2a) 

,2 N „i »i n 2t/ N „i n 

^«+i-yE / ^^'E'^^''(c"(^)-cK.'))-^ E / ds'l[6('^\c'^is)-c^,is')). (2.2b) 

j.j' — l ^ a— j.j' — l ^ a— 

Here, (• • •)^]^ denotes a thermal average taken with respect to the Hamiltonian for n + 1 replicas of the noninter- 
acting, uncrosslinked system of macromolecules. Moreover, Ti^^^ effective pure Hamiltonian accounting for the 
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interactions amongst the macromolecules and the effects of crosshnking, the latter generating interactions between 
the replicas. The parameter measures the strength of the excluded-volume interaction; the parameter /x^ measures 
the density of the constraints and serves as the control parameter for the vulcanization transition. As a result of there 
being random constraints rather than interactions, the coupling between the replicas takes the form of product over 
all replicas rather than, say, a pairwise sum. As usual, the disorder-averaged free energy is proportional to [InZ], 
which we obtain via the replica technique as lim„_^o In [^"] ■ Let us mention, in passing, the symmetry content 
of this replica theory: Ti-fi^i is invariant under arbitrary independent translations and rotations of the replicas as well 
as their arbitrary permutation. 

The natural collective coordinates for the vulcanization transition are 

^(^) = ]V^/ e'^P*^ • (2-3) 

3 = 1 "'0 



which emerge upon introducing Fourier representations of the two types of delta function in Eq. ( ^.2b ), as discussed in 



detail in see Sec. 5.1 of Ref. ^ . (Such collective coordinates were first introduced in the context of crosslinked macro- 
molecular melts by Ball and Edwards |^^.) We use the symbol k to denote the replicated wavevector {k°, k^, . . . , k"}, 
and define the extended scalar product fc • c by k" • c" -I- k^ • -I- • • • -I- k" • c". The collective coordinates Q{k) are the 
microscopic prototype of the order parameter (pH]), the latter being related to Q{k) via hm„^o {Q{k))^^j^, where 

, .p _ (■■•exp(-HP+i))^^^ 
(exp(-H,^+i)) 



C. Replica field theory for vulcanized macromolecular systems 

As discussed in detail in Sec. 5.3 of Ref. [Q, one can put the partition function into a form of a field theory by 
applying a Hubbard- Stratonovich transformation to the collective coordinates Q{k); we denote the corresponding 
auxiliary order-parameter field by ^(k). At this stage one encounters a vital issue, viz., that it is essential to draw 
the distinction between examples of Q{k) and ^l{k) that belong to the one-replica sector (IRS) and those that belong 
to the higher-replica sector (HRS). The distinction lies in the value of k: for a replicated wavevector fc, if there is 
exactly one replica for which the corresponding d- vector k" is nonzero [e.g. k — (0, . . . , 0, k" ^ 0, 0, . . . , 0)] then we 
say that k lies in the one-replica-sector (fc e IRS) and that the corresponding Q{k) and il{k) are IRS quantities. On 
the other hand, if there is more than one replica for which the corresponding components of k are nonzero then we 
say that k lies in the higher-replica-sector (k £ HRS) and that the corresponding Q{k) and U{k) are HRS quantities. 
For example, if fc = (0, . . . , k" 7^ 0, . . . , k*^ 7^ 0, . . . , 0) then k lies in the HRS. (More specifically, in this example k 
lies in the two-replica sector of the HRS.) The importance of this distinction between the IRS and the HRS lies in 



the fact, evident from the order parameter (2T), that the vulcanization transition is detected by fields residing in the 
HRS, whereas the IRS fields measure the local monomer density, and neither exhibit critical fluctuations near the 
vulcanization transition nor acquire a nonzero expectation value in the amorphous solid state. 

Bearing in mind this distinction between the IRS and HRS fields, the aforementioned Hubbard-Stratonovich trans- 
formation leads to the following field-theoretic representation of the disordered-averaged replicated partition function: 

[Z"] cx J Vln j V^n exp ( - ndNJ'n{{^°'{k),n{k)})y (2.5) 

where 17" (k) [which represents n{k) when fc = (0, . . . , 0, k" = k 7^ 0, 0, . . . , 0)] is a IRS field, n{k) is a HRS field, and 
the explicit expressions for the resulting effective Hamiltonian J-n and functional integration measures are given by 
Eqs. (5.12) and (5.9) of Ref. [^. In this formulation of the statistical mechanics of RCMSs, one can readily establish 
exact relationships connecting average values and correlators of Q{k) with those of Q{k) [^. (Such relationships 
between expectation values involving microscopic variables and auxiliary fields are common in the setting of field 
theories derived via Hubbard-Stratonovich transformations [p6|.) For example, for wavevectors lying in the HRS one 
has 
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2jY "fe+fc',0 ' 



(2.6a) 
(2.6b) 



where {• ■ denotes an average over the field the ory (2T5| ), a nd th e subscript c mdicates that the correlators are 

connected. Relationships such as those given in Eqs. ( 2.6a ) and ( 2.6b ) allow one to relate order-parameter correlators 
to correlators of the field theory. 



D. Minimal model for the vulcanization transition 



The exact field-theoretic representation of RCMSs discussed in the previous section serves as motivation for a 
minimal model capable of describing the universal aspects of the vulcanization transition inasmuch as it indicates 
the appropriate order parameter and symmetry content. In the spirit of the standard Landau approach, one can 
determine the form of the minimal model by invoking symmetry arguments along with three further assumptions: 
(i) that fluctuations representing real-space variations in the local density of the constituents are free-energetically 
very costly, and should therefore be either suppressed energetically or, equivalently (as far as our present aims are 
concerned), prevented via a kinematic constraint; (ii) that we need only consider order-parameter configurations 
representing physical situations in which the fraction of constituents localized is at most small; and (iii) that the field 
components responsible for the incipient instability of the liquid phase are those with long wavelengths. Provided these 
assumptions hold, one may: (i) expand the effective Hamiltonian in powers of the order parameter; and (ii) expand 
the coefficient functions in powers of wavevectors. One retains terms only to the order necessary for a description of 
both sides of the transition. (When we go beyond mean-field theory, below, KG arguments will justify our omission of 
all other symmetry-allowed terms on the grounds that they are irrelevant at the fixed-points of interest.) This scheme 
leads to the following minimal model ||l^,|2^, which takes the form of a cubic field theory involving a HRS field il(fc) 
that lives on (n + l)-fold replicated d-dimensional space: 

[Z"] cx y P^rjcxp(-5„), (2.7a) 
S^{m)^N [-aT+^\k\^)\n{k)\^^Ng ^ ^ik) n{k) S^^^f^^^j^^^^ , (2.7b) 

fceHRS fci,fc2,fe3eHRS 

where r is the reduced control-parameter measuring the crosslink density. This model was introduced in Ref. |ic| ] 
as a Landau theory of the vulcanization transition, where it was shown to yield a rich description of the amorphous 



solid state, even at the saddle-point level, which we briefiy summarize in Sec. [II (along with the results of various 
semi- microscopic approaches). Although the semi- microscopic derivation of 5„ containsn-dependent coefficients a„, 
bn and it is admissible for us to keep only the n — > limit of these coefficients (i.e. a, b and g) at the outset 
because iS„ is already proportional to n for pertinent field-configurations . 

We wish to emphasize the point th at th is minimal model does not contain fields outside the HRS. For example. 



in the cubic interaction term in Eq. (2.7b), the wavevectors in the summations are constrained to lie in the HRS. 
This (linear) constraint on the field embodies the notion that inter-particle interactions give a "mass" in the IRS 
(i.e. produce a free-energy penalty for density inhomogeneities) that remains nonzero at the vulcanization transition. 
From the standpoint of symmetry, this constraint has the effect of ensuring that the only symmetry of the theory 
(associated with the mixing of the replicas) is the permutation symmetry S„+i. Without it, the model would have 
the larger symmetry, 0((n-|- l)d), of rotations that mix the (Cartesian components of the) replicas; see the term 
associated with the inter-replica coupling arising from the disorder-averaging of the replicated crosslinking constraints 



in Eq. ( p. 2b ). In addition to permutation symmetry, the model has the symmetry of independent translations and 



rotations of each replica. The restriction to the HRS (or, equivalently, the energetic suppression of the IRS) is vital: 
it entirely changes the content of the theory. Without it, one would be led to completely erroneous results for both 
the mean-field picture of the amorphous solid state and, as we shall see, the critical properties of the vulcanization 
transition. 



For use in Sec. V A, when we come to examine the physical implications of the Ginzburg criterion, we list values of 
the coefficients in the action derived for the case of RCMSs (up to inessential factors of the crosslink density control 
parameter fj,^): 

T^{^i'- (2.8a) 

a =1/2, (2.8b) 
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b = M/6d, (2.8c) 
g = 1/6. (2.8d) 

Here, /i^ is the mean- field critical value of /i^, i is the arclength of each macromolecule, and i is the persistence length 
of the macromolecules. 

III. VULCANIZATION TRANSITION IN MEAN-FIELD THEORY: BRIEF SUMMARY OF RESULTS 
A. Mean-field order parameter: Liquid and amorplious solid states 



Mean- field investigations of RCMSs and related systems |3|,[4|,|10|-|12| have shown that: (i) There is a continuous 
phase transition between a liquid and an amorphous solid state as a function of the density of the crosslinks (or 
other random constraints). This transition is contained within the HRS. Both the liquid and the amorphous solid 
states have uniform densities, and therefore the order parameter is zero in the IRS on both sides of the transition, 
(ii) In the solid state, translational invariance is spontaneously broken at the microscopic level, inasmuch as a nonzero 
fraction of the particles have become localized in space. However, owing to the randomness of the localization, this 



symmetry-breaking is hidden. [Hence the need for a subtle order parameter (2.1).] In the language of replicas, the 
symmetries of independent translations and rotations of the replicas are spontaneously broken, and all that remains 
are the symmetries of common translations and rotations (corresponding to the macroscopic homogeneity and isotropy 
of the amorphous solid state). The permutation symmetry amongst the n -I- 1 replicas appears to remain intact at 
the transition, (iii) The stationarity condition for the order parameter can be solved exactly. In the context of the 
minimal model, in the liquid state one finds fl{k) — 0; in the solid state the order parameter takes the form 



n{k) = {2aT/3g)S^^l w(y^afcV6r), (3.1a) 
w{k)= / deT:{e)e~'^^'^\ (3.1b) 



where k = X)q=o^"- function 7r(6') is a universal function, in the sense that it does not depend on the model- 
specific coefficients a, b and g: it is normalized to unity and satisfies a certain nonlinear integro-differential equation; 
see Refs. ||^,^,|l0| . From the physical perspective, uj{k) encodes the distribution of localization lengths of the localized 
monomers and the Kronecker delta factor 5^^^^ exhibits the macroscopic translational invariance of the random solid 



state. By passing to the k Q limit in Eq. (B.la) one learns that the fraction of localized monomers q (i.e. the gel 
fraction) is given by 

f 0, liquid state; ,„ 

^ " \ (2a/3.g) r'', solid state; > 

with the exponent /3 being given by the mean-field value of unity. It has recently been demonstrated that the mean-field 
state summarized here is locally stable [28| . (We note, in passing, that no spontaneously replica-symmetry-breaking 
solutions of the order-parameter stationary condition have been found, to date.) 

B. Gaussian correlator: Liquid and critical states 

The incipient amorphous solidification, as the vulcanization transition is approached from the liquid side, is marked 
by strong order-parameter fluctuations, which are diagnosed via the correlator G{k) defined through 

^ ^(^'))r+i,c • (3-3) 



The unusual factor of 1/iV is due to our choice of the normalization of Q{k) in Eq. ( |2.3|). S ection IV, below, is 
dedicated to explaining the physical content of this correlator and precisely how, via Eq. ( 2.6bD , it is able to detect 



incipient random solidification. The value of the correlator in the mean-field approximation follows from the quadratic 



terms in Eq.( 2.7b) and is given by 
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G(fc) « G,{k) ^ ^ ^ , (3.4) 

which below wiU play the role of the bare propagator. Notice that Go{k) obeys the homogeneity relation 

GCk,T)^\k\-'+^g{\k\\rr), (3.5) 

in which g{x) ^ x^^"^ for x +0 and approaches a constant value for large x. Moreover, the exponents take on the 
mean-field values rj — Q, u — 1/2 and 7 = v{2 — 77) = 1, this last relationship guaranteeing that the susceptibility 
limj,^g G(A;, t) diverges as \t\^'^ ■ 



IV. ORDER-PARAMETER CORRELATOR AND SUSCEPTIBILITY, AND THEIR PHYSICAL 

SIGNIFICANCE 



Let us now consider the order-parameter correlator and the associated susceptibility from the perspective of incipient 
random localization |2^]. In the simpler context of, e.g., the ferromagnetic Ising transition the two-point spin-spin 
correlator quantifies the idea that the externally-imposed alignment of a particular spin would induce appreciable 
alignment of most spins within roughly one correlation length of that spin, this distance growing as the transition is 
approached from the paramagnetic state. How are these ideas borne out in the context of the vulcanization transition? 
Imagine approaching the transition from the liquid side: then the incipient order involves random localization and so, 
by analogy with the Ising case, the appropriate correlator is the one that addresses the question: Suppose a monomer 
is localized to within a region of some size by an external agent: Over what region are other monomers likely to 
respond by becoming localized, and how localized will they be? We can also consider the order-parameter correlator 
and the associated susceptibility from the perspective of the formation of (mobile, thermally fluctuating) assemblages 
of macromolecules, which we refer to as clusters: How do they diagnose the development of larger and larger clusters 
of connected macromolecules, as the crosslink density is increased towards the vulcanization transition? 

Bearing these remarks in mind, we now examine in detail the physical interpretation of the order-parameter cor- 
relator {Q{k) Q(~fc))„_|_]^ which, as we shall see, captures the physics of incipient localization and cluster formation. 
To see this, consider the construction 



Ct(r-r') 



V 

N 



J2 ds ds'(5W(r-c,(.))<5W(r'-c,,(.'))) 

-.ji^lJo Jo 

X ( exp —it ■ {cj (s) — r) exp it ■ {cj'(^s') ^ r') 



(4.1a) 



which, in addition to depending on the separation r — r', depends on the "probe" wavevector t. The first expectation 
value in this construction accounts for the likelihood that monomers {j,s) and {j',s') will respectively be found 
around r and r'; the second describes the correlation between the respective fluctuations of monomer (j, s) about r 
and monomer (j',s') about r'. 

Now, the quantity Ct(r — r') is closely related to an HRS correlator involving the semi- microscopic order parameter 
Q{k). To see this we introduce Fourier representations of the two delta functions and invoke translational invariance, 
thus establishing that [^o) 



,i(k+t)-(r-r') 



V 

N 



Ar2 



_ g.k.(r-r') ^^^^ k - t, t, 0, . . . , 0)* Q(0, k - t, t, 0, . . . , 0))„^^. 



(4.2) 
(4.3) 



Having seen that Ct(r — r') is closely related to an HRS correlator involving Q(}i) [which can be computed via the 
fi field theory], we now explain in more detail how Ct(r — r') detects the spatial extent of relative localization. First, 
let us dispense with the case of t = 0. In this case Ct(r — r') is simply {V/N times) the real-space density-density 
correlation function and, as such, is not of central relevance at the amorphous solidification transition. Next, let us 
consider the small-t limit of Ct(r — r'). This quantity addresses the question: If a monomer at r is localized "by hand" 
what is the likehood that a monomer at r' responds by being localized at all, no matter how weakly. It is analogous 
to the correlation function defined in percolation theory that addresses the connectedness of clusters 1 21 . 
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To substantiate the claim made in the previous paragraph we examine the contribution from each pair of monomers 
to the quantity Ct(r — r'). Let us start from the simplest situation, in which no crosslinks have been imposed. We 
assume that t is small (i.e. V^-^^^ '3> |t|^^ ^ i?g, where Rg is the radius of gyration for a single macromolecule) 
and that the macromolecular system has only short-range interactions. For each term in the double summation over 
monomers there are two cases to consider, depending on whether or note the pair of monomers are on the same 
macromolecule. For a generic pair of monomers that are on the same macromolecule (i.e. j = j'), we expect that 
(expit • (cj(s) - Cj{s'))) - 1, and that (for |r - r'| < Rg) (S'^'^'>{r - Cj{s)) ,5W(r' - Cj{s'))) ~ V'^ R-<^. Then the 
total contribution to Ct(r — r') coming from pairs of monomers on the same macromolecule is of order {N/V)'^R~'^. 
On the other hand, for a generic pair of monomers that are on different macromolecules (i.e. j 7^ j'), we expect that 
(expit • (cj(s) - Cj,{s'))) ~ and that (^(''^(r - Cj(s)) 5('^)(r' - Cj,{s'))) ~ V-"^. Therefore the total contribution 
to Ct(r — r') coming from pairs of monomers on different macromolecules is of order {N/V)^V~^ . Thus, we find that 
the intrachain (i.e. j — j') contribution to Ct(r — r') dominates over the interchain (i.e. j 7^ j') contribution in the 
thermodynamic limit. 

Moving on to the physically relevant case, in which crosslinks have been introduced so as to form clusters of 
macromolecules, we see that what were the intrachain and interchain contributions become intracluster and intercluster 
contributions. With the appropriate (slight) changes, the previous analysis holds, which indicates that the intracluster 
contribution dominates Ct(r — r') in the thermodynamic limit. In other words, in the small-t limit a pair of monomers 
located at r and r' contribute unity to Ct(r — r') if they are on the same cluster and zero otherwise. This view allows 
us to identify the small-t limit of Ct(r — r') with the pair-connectedness function defined in (the on-lattice version of) 
percolation theory |^ . 

What about Ct(r — r') in the case of general t? In this case it addresses the question: If a monomer near r is 
localized on the scale (or more strongly), how likely is a monomer near r' to be localized on the same scale (or 
more strongly)? This additional domain of physical issues associated with the strength of localization results from 
the effects of thermal fluctuations, and is present in the vulcanization picture but not the percolation one. 

Let us illustrate the significance of Ct(r — r') by computi ng it in the setting of the Gaussian approximation to 



the liquid state in three dimensions. To do this, we use Eq. ( 2.6b ) to express Ct(r — r') in terms of the (Gaussian 

)n+l,c' 



approximation to the) correlator (f2(fc) il(fc'))^^]^ ^, which has the Ornstein-Zernicke form given in Eq. ( |3.4| ). Thus, 



we arrive at the real-space Yukawa form 



|Ct(r-rO|oc^^^Pl±l4f^, (4.4a) 
I I |r — r I 

where the correlation length C is defined by = — 2ar. Hence, we see the appearance of a probe-wavelength- 
dependent correlation length Coff(i)- The physical interpretation is as follows: in the t ^ limit, Ct(r — r') is testing 
for relative localization, regardless of the strength of that localization and, consequently, the range of the correlator 
diverges at the vulcanization transition. This reflects the incipience of an infinite cluster, due to which very distant 
macromolecules can be relatively localized. By contrast, for generic t it is relative localization on a scale t^^ (or 
smaller) that is being tested for. At sufficiently large separations, even if a pair of macromolecules are relatively 
localized, this relative localization is so weak that the pair does not contribute to Ct(r — r'). This picture is reflected 
by the fact that Ceff(0 remains finite at the transition. 

Given that we have identified a correlator that is becoming long-ranged at the transition, it is natural to seek an 
associated divergent susceptibility 9t. To do this, we integrate Ct(r — r') over space and obtain 

et= / ^^^^-^Ct(r-r') = iVlim (Q(0,t,-t,0,...,0)*g(0,t,-t,0,...,0))|;^,. (4.5) 
Passing to the t ^ limit, we have 

limOt- (-r)--^, (4.6) 
t^o 

where the final asymptotic equality is obtained from a computation of the (field-theoretic) corr elator {Vlik) Vlik')) ^ 

[see Eq. ( ^.6b| )]. This quantity is measure of the spatial extent over which pairs of monomers are relatively localized, 
no matter how weakly, and thus diverges at the vulcanization transition. At the Gaussian level of approximation. 



Eq. (3.4), this susceptibility diverges with the classical exponent 7=1. By contrast, for generic t the susceptibility 
Ot remains finite at the transition, even though an infinite cluster is emerging, due to the suppression of contributions 
to Ot from pairs of monomers whose relative localization is sufficiently weak (i.e. those that lead to the divergence in 
the small-t limit). 
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V. VULCANIZATION TRANSITION BEYOND MEAN-FIELD THEORY 



A. Ginzburg criterion for the vulcanization transition 



To begin the process of analyzing the vulcanization transition beyond the mean-field (i.e. tree) level, we estimate 
the width St of reduced constraint-densities r within which the effects of order-parameter fluctuations about the 
saddle-point value cannot be treated as weak, i.e., we construct the Ginzburg criterion. To do this, we follow the 
conventional strategy (see, e.g., Ref. [3l| ]) of computing a loop expansion for the 2-point vertex function to one-loop 
order and examining its low-wavevector limit (i.e. the inverse susceptibility). Note that in the present setting the loop 
expansion amounts to an expansion in the inverse monomer density. 

A A 



ki - k 



A ^\ A 

k X \ k 




FIG. 1. One-loop correction to the 2-point vertex function. Full lines indicate bare HRS correlators; dashed lines indicate 
amputated external bare HRS correlators. 



Our starting point is the minimal model, Eq. ( p. 7b ), for which the bare correlator is given by Eq. (3.4). Then the 



one-loop correction to the 2-point vertex function comes from the diagram shown in Fig. |^, which is calculated in 
App. 0. By choosing k e 3+RS ( i.e. in the HRS but not in the two-replica sector) |p2| we obtain for the inverse 
susceptibility the result 

(«xr'^-2ar+184/p^-^, (5^1) 

in which a large wavevector cut-off at |fc| = A is implied. The (one-loop) shifted critical point Tc marks the vanishing 
of x~^, i-e., solves 

2V f d'^P 



= -2are + 18<7^— / — ^ , . (5.2) 

Now, in mean-field theory the transition occurs a t t — 0, with positive (resp. negative) values corresponding to the 
amorphous solid (resp. liquid) states. From Eq. ( |5.2| ) we see that that inclusion of fluctuations enlarges the region 
of crosslink densities in which liquid state is stable, as one would expect on general physical grounds. However, it is 
worth noting, in passing, that without the exclusion of the one-replica sector the converse would occur (i.e. fluctuations 
would enlarge the region of stability of the amorphous solid state). By subtracting Eq. ( |5.2| ) from Eq. (^) in the 
standard way, replacing Tc by its mean-field value (of zero) in the loop correction, and rescaling the integration variable 
according to — —2aTk'^, we arrive at 

(Nx)-' = -2a{T - r,) (l - l8g'{V/N)b-^/^{~2aT)^^-'^^/\Ja) , (5.3) 

where i s a dimensionless number dependent on d (and weakly on A, at least in regime of interest, i.e., d below 6). 
Equation ( |5.3| ) shows that for d < 6 a fluctuation dominated-regime is inevitable for sufhcient small r, and hence that 
the upper critical dimension for the vulcanization transiti on is six, in agreement with naive power-counting arguments 
applied to the n —> limit of the cubic field theory, Eq. ( p. 7b ). The Ginzburg criterion amounts to determining the 



departure of r from its critical value such that in Eq. ( ^.3| ) the one-loop correction is comparable in magnitude to the 
mean- field level result. 
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To determine the physical content of the Ginzburg criterion, we i nvoke t he va lues of the coefficients of the minimal 
model appropriate for the semi-microscopic model of RCMSs, Eqs. ( 2.8a )-( 2.8d ), and we exchange the macromolecule 
density N/V for the volume fraction = {N/V){L/£)£'^. Thus we arrive at the following form of the Ginzburg 
criterion: for d < 6, fluctuations cannot be neglected for values of r satisfying 



< 



[L/ey 



(5.4) 



from which we sec that the fluctuation-dominated regime is narrower for longer macromolecules and higher densities 
(for 2 < d < 6). Such dependence on the degree of polymerization L/i is precisely that argued for long ago by 
de Gennes on the basis of a percolation-theory picture |Q . 

Besides the fields and vertices featuring in the minimal model, there are other symmetry-allowed fields and vertices 
that are generated by the semi-microscopic theory of RCMSs. Examples are provided by the IRS field, which describes 
density fluctuations, along with vertices of cubic, quartic or higher-order that couple the IRS field to the HRS field. 
In App. ^ we investigate the effect of these fields and vertices, which are omitted from the minimal model, and 
show: (i) that the inclusion of their effects (at the one-loop level) does not change the Ginzburg criterion derived in 
the present section; and (ii) that the HRS critical fluctuations do not provide any singular contributions to the IRS 
density-density correlation function (at least to one-loop order). 

B. Renormalization-group procedure and its subtleties 

We now describe the RG procedure that we are using, a schematic depiction of which is given in Fig. |. The mam 
thrust of our approach is the standard "momentum-shell" RG, via which we aim to determine how the parameters of 
the theory, t and g, flow under the two RG steps of coarse graining and rescaling. However, in the present context 
there are some significant subtleties owing to the need to constrain the fields to lie in the HRS. 

In the coarse-graining step, we integrate out the rapidly-varying components of fl(k) (i.e. those corresponding to 
wavevectors satisfying A/b < \k\ < A). Here, the constraint that only the HRS field is a critical field demands that 
one treat the HRS and the IRS distinctly. We handle this by working with a large but finite (replicated) system 
contained in a hyper-cubic box of volume 1/"+^ on which periodic boundary conditions are applied. As a consequence, 
the wavevectors are "quantized" and therefore we can directly make the appropriate subtractions associated with the 
removal of the zero- and one-replica sectors. Having made the necessary subtractions, we compute the various Feynman 
diagrams (for the construction of the Ginzburg criterion and the the coarse-graining step of the RG) by passing to 
the continuous wavevector limit (so that wavevector summations become integrations) . 

The replica technique has the following curious feature. In the infinite-volume limit the different sectors are spaces 
of different dimensionalities, and thus the contributions from the lower replica sectors appear to be sets of measure 
zero relative to the contributions from the HRS. However, in the replica limit, the contributions from different sectors 
are comparable and, hence, the lower sectors cannot be neglected. 





h- 



A/b 

— \ — 



o I I I I I I I I I I I I I I I i f i I I I I I h ] (a) 



o I I I I I I I I I I I I I I I I (b) 

\\\ 

o I I I I I I I I I I I I I I I I (c) 



o I I I I I I I I I I I I I I I I I I I I I I I I (d) 
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FIG. 2. Schematic one-dimensional depiction of the basic steps of the RG procedure (the field variables are defined only at 
the hash marks denoting the quantized wavevectors) : from (a) to (b) integrate out the fields at the quantized wavevectors k in 
the "momentum shell" (shaded) ; from (b) to (c) rescale lengths to restore the wavevector cut-off, and rescale the field to restore 
the gradient term; from (c) to (d) restore the density of the degrees of freedom. (In practice, we employ a momentum shell of 
infinitesimal width.) 



The coarse-graining step is followed by the rescaling step, in which the aim is to return the theory to its original 
form. The field- and length-rescaling aspects of this step (to recover the original wavevector cut-off and form of 
the gradient term) are standard, but there is a subtlety associated with the fact that the original theory is defined 
on a finite volume (in order that the wavevectors be quantized and the various replica sectors thereby be readily 
identifiable). This subtlety is that upon coarse-graining and rescaling one arrives at a theory that is almost of the 
original form, but is defined on a coarser lattice of quantized wavevectors associated with the reduced (real-space) 
volume. If we wish to return the theory to its truly original form, we are required to increase the density of the 
coarsened wavevector lattice. To accomplish this, we choose to make use of the extension to (n -|- l)d dimensions of 
the following one-dimensional relation, exact in the thermodynamic (i.e. large real-space size B) limit: 

ke{2TTnb/B} ke{2TTn/B} 



One way to understand this is to regard the two sides of Eq. (5.5) as providing different discrete approximations 



to the same continuous-wavevector (i.e. infinite-volume) limit. Thus, we expect the difference bet ween them to be 



unimportant in the thermodynamic limit. Another way is to regard the right hand side of Eq. (5.5) as pertaining to 
a system with a larger number of degrees of freedom than the left hand side, but that the factor appropriately 
diminishes the weight of each degree of freedom. It would be equally satisfactory if we chose, in our RG scheme, not 



to restore the wavevector lattice spacing, which would amount to our using the left-hand-side of Eq. (|5.5|) 



C. Expansion around six dimensions 

In the previous two subsections we have established that the upper critical dimension for the vulcanization transition 
is six, and we have described an RG procedure capable of elucidating certain universal features of the transition. We 
now examine the RG flow equations near the upper critical dimension that emerge from this procedure, and discuss 
the resulting fixed-point structure and universal critical exponents. To streamline the presentation we have relegated 
the technical details of the derivation of the flow equations to App. O 



1. Flow equations 



As with the mean-field theory and the Ginzburg criterion, our starting point is the replicated cubic field theory, 
Eq. ( 2.7b| ). By suitably redefining the scales of n{k) and k we can absorb the coefficients a and b, hence arriving at 
the Landau- Wilson effective Hamiltonian 

S„{{n})^N J2 (-r+l\k\')\n{k)\' -Ng ^(k,) n{k) Sf,^^^,^^,,^^^ , (5.6) 



fceHRS /ci,fe2,fe36HRS 

in which all wavevector summations are cut off beyond replicated wavevectors of large magnitude A, from which we 
can read off the bare correlator 

Go(fc) = — ■ (5.7) 
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(a) 



(b) 



FIG. 3. Contributing one-loop diagrams. Full lines indicate bare HRS correlators for short-wavelength fields (i.e. fields lying 
in the momentum shell); wavy lines indicate long- wavelength fields. 

We shall be working to one-loop order and, correspondingly, the diagrams that contribute to the renormahzation 
of the parameters of the Landau- Wilson effective Hamiltonian are those depicted in Figs. || (a) and (b). The resulting 
flow equations are 



dg/d\nb = g{e/2 ~ C^g^ - ^C,g^ + 0{Tg\eg\ g^)) , 
dzld\nb=]^{d + 2-Cig^)+0{Tg\eg\g^), 



(5.8a) 
(5.8b) 

(5.8c) 



where £ = 6 — d, 6 is the length-rescaling factor, z is the field-rescaling factor, and the (constant) coefficients in the 
flow equations are given by 



(Oo, Oq, Oi, O3) — 



(9A^36,-6,72), 



TV (27r)6 

in which S% is the surface area of a 6-dimensional sphere of unit radius. 



(5.9) 



2. Fixed-point analysis and its consequences 



We proceed in the standard way by first finding the fixed points (t* , g^) of the flow equations, at which o?(t, g)/dhib 
(0,0). We linearize the flow equations about each of the resulting fixed points. 



din 6 V 5^5^ 



2 - (C^ + 



~2Co 5* 



(5.10) 



where we have dropped higher-order corrections. We then establish the RG eigenva lues at each fixed point by finding 
the eigenvalues of the linearized RG transformation matrices. Finally, we solve Eq. ( 3.1C| ) to obtain the flow near each 
fixed point. 

For e both negative and positive (i.e. for d both above and below six) we find a Gaussian fixed point (GFP): 
(T,, g,) = (0,0). Solving Eq. (5.10) about this fixed point gives the flow 



g{b) ~ [g{l)by^ 



(5.11) 
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with the RG eigenvalues yi and 1 /2 respectively given by j/r = 2 and Ug = e/2. 

As one can see from Eq. ( ^.11 ), above six dimensions the GFP is unstable in the r direction and stable in the g 



direction. However, below six dimensions the GFP also becomes unstable in the g direction, and a new fixed point- 
the Wilson-Fisher fixed point (WFFP)— emerges, located at {n,gl) = ((AV28), (l/126)((27r)V5'6)(y/A^)-^) e. (Let 
us mention, in passing, that if we had not correctly implemented the constraint that wavevector summations exclude 
contributions for the IRS then the structure of the flow equati ons w ould have been utterly different; e.g., the WFFP 
would have occurred at a complex value of g.) By solving Eq. ( |5.10| ) for the WFFP we find the flow 

where yl= {3/VT4){{V/N){Se/{2n)^)y^'^ {A'^e^/'^) and the RG eigenvalues are given by yi = 2- (5£/21) and = 

We now proceed to obtain the critical exponents for physical quantities from the RG eigenvalues at each fixed point. 
The homogeneity relation for the correlator G{k), following from a standard RG analysis |^^, reads 

G(fc, r) = 02 ^-d Qi^fj^^ ^y^^y (5 ^^3^ 

We eliminate h by choosing h\k\ = 1; then comparison with Eq. (3.5) leads to ly — 1/yi and 77 = Cig^. Thus, for the 



GFP we have 

iy-^ = 2, ri = 0, (5.14) 

and for the WFFP we have, to first order in e, 

ly-^ =2- {5e/2l), rj = -e/2l. (5.15) 

Both above and below six dimensions, the critical exponents ly and 77 (and (3, to be discussed below) are identical to 
those governing analogous quantities in percolation theory (at least to first order in e), as computed via the Potts 
field theory l^^l- We discuss the significance of this result and the relationship between the present approach and 
percolation/gelation-based approaches in Sec. [v| . 

We have focused on the cubic interaction in the vulcanization field theory. There are, of course, additional symmetry- 
allowed interactions, such as the quartic interaction. Near to six dimensions, however, the fact that such interactions 
are irrelevant at the GFP can be shown by naive power-counting arguments, which hold in the replica limit (and 
remain uncompromised at the WFFP, owing to its proximity to the GFP). 



D. Scaling for gel fraction and wavevector-dependent order parameter 

In order to relate properties of the amorphous solid state to those computed in the liquid and critical states, we now 



follow the standard scahng analysis. To do this, we add to the minimal model, Eq. (5^), a source field that couples 
linearly to the order parameter: —Nj^kews ^(^) U{—k). We assume that U contains only long wavelength compo- 
nents, so that it does not couple to any field featuring in any momentum-shell integrations. Then the renormalization 
of U{—k) comes only from the rescalings of k and ri(fc), and thus we have 

U'{k')^ zb-'^U(k). (5.16) 

To obtain the exponent /3, which describes the scaling of the gel fraction the conventional method prescribes the 
application of a uniform source field. In the present theory, the (zero replica sector) field variable ri(0), which would 
couple to such a uniform source, is excluded, and instead we choose U{k) ~ f^^k+ko 6' '^l^'^re fco lives in the HRS but 
is otherwise arbitrarily small. (This prescription is consistent with the notion that the gel fraction follows from the 
long- wavelength limit of the order parameter, the limit being taken via wavevectors in the HRS.) Hence we arrive the 
recursion relation for h: 

h'^zh^by^- yh = id + 2-ri)/2. (5.17) 

As we are already in possession of 77 at the GFP and the WFFP, we thus arrive at the scaling dimension y^ of the 
source field h. 

Having obtained yn, we now use it, together with y^, yg and the singular part of the free energy density /, to 
determine (3, in the following way. According to homogeneity, / has the form 
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/(r, g, h) = b-^firby- , gby^ , hby^ ) . (5.18) 

By taking the derivative with respect to h so as to form the order-parameter equation of state, choosing h = 0, and 
passing to the small fco limit, one finds the following scaling behavior of the gel fraction: 

g(T,5,0)- lim - 6"'^+^'' M(T6^-,g5^«,0) = T(''-^'-)/^-Af(l,.gr-^<'/^-,0). (5.19) 

feo->6 

Let us first consider the regime d > 6, for which the appropriate fixed point is the GFP and, therefore one expects 
the exponents to take on their classi cal va lues. Now, as one can see from the mean-field value for the order parameter 
n (and thus the gel fraction q), Eq. ( j^.la ), both of which are proportional to g~^, the cubic interaction is dangerously 



irrelevant at the GFP, and thus one has 



Hence, near the GFP one has 



M(l,5,0)-^, for g->+0. (5.20) 



Q(T,.g,0)~T'3, for g^+0, (5.21a) 

I J d+2 I 6-d 

/3 = + ^ = + ~ = 1, (5.21b) 

Vt Vt 2 



which is precisely the mean- field value of the exponent (3 given in Sec. [Ill A . 

Now let us turn to the regime d < 6, for which the exponents are nonclassical. The appropriate fixed point is now 
the WFFP, at which the cubic interaction is irrelevant but not dangerously so. Thus, in this regime one has the 
standard scaling relation 

/3 = ^^ = l-(£/7), (5.22) 

Vt 

where the second equality holds only to order e. 

In fact, under the (not unreasonable) assumption that there is only one characteristic length-scale in the ordered 
state (i.e. that the fluctuation correlation length does not provide a length-scale independent from the localization 
length-scale), we can go beyond the establishing of the scaling of the gel fraction (i.e. the long- wavelength limit of the 
order parameter) and propose a more general scaling hypothesis, which incorporates the scaling of the (singular part 
of the) wavevector-dependent order parameter p7[ |. This takes the form of the scaling hypothesis: 

{n{k)) o^t^w^Pt^''). (5.23) 

The quantity r"^, which plays the role of the fluctuation correlation length in the liquid state, is here seen to play the 
role of the characteristic scale for the localization lengths in the ordered state. Presumably, it also governs the scale 
over which (amplitude-type) fluctuations are correlated in the solid state. Let us note that the mean-field result for 
the order parameter not only obeys this scaling relation (with j3 = 2v = 1) but also provides an explicit form for the 
function w. 

VI. CONCLUDING REMARKS: CONNECTIONS WITH OTHER APPROACHES AND THE ROLE OF 

THERMAL FLUCTUATIONS 

Having constructed an RG theory for the liquid and critical states of vulcanized matter, we now examine the 
results of this RG theory and discuss the relations hip be tween these results and the results of other approaches to the 



vulcanization transition. As we have seen in Sees. V C , via an expansion around six spatial dimensions our minimal 
model for the vulcanization transition yields values for certain critical exponents that characterize the behavior of 
the system near to and at the transition. These exponents turn out to be numerically equal to those characterizing 
physically analogous quantities in percolation theory, at least to first order in the departure e from six dimensions. 
We have not proven that the equality between exponents holds beyond first order in e, although there are hints in 
the structure of the theory suggesting that it does. 

This equality between exponents seems reasonable in view of the intimate relationship between percolation theory 
and the connectivity of the system of crosslinked macromolecules, this connectivity pertaining to the statistics of 
systems formed according to the Deam-Edwards distribution of quenched randomness (and hence to the statistical 
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mechanics of the uncrosshnked macromolecular liquid) [Q. Indeed, a connection between the percolation and vul- 
canization transitions already shows up at the level of mean-field theory: the dependence of the gel fraction q on the 
crosslink-density control parameter fM^ obtained via the semi- microscopic approach (in the case of RCMSs), viz., that 
q obeys 

1 - g = cxp (-/i^g) , (6.1) 

is identical to the mean-field-percolation dependence of the fraction of sites participating in the infinite cluster, 
obtained by Erdos and Renyi in their work on random graphs [ |39| , this identity holding not just near the transition, 
where the dependence of q on /i^ — linear, but for all crosslink densities. Moreover, the mean-field result emerging 
from the minimal model of the vulcanization transition yields this linear dependence (but cannot, of course, be applied 
beyond the transition regime) . The relevance of percolation theory to the vulcanization transition also manifests itself 



beyond the mean-field level in the physical meaning of the order-parameter correlator, as we have discussed in Sec. IV. 
This connection has long been realized, and supports the use of percolative approaches as models of certain aspects 
of the vulcanization transition |p^|-^o|] . 

These percolative approaches include direct applications of percolation theory JT6|-p^ , mentioned in the preceding 
paragraph, as well as the approach given by Lubensky and Isaacson ^o| . The latter approach extends the connection 
between the statistics of linear macromolecules and the zero-component limit of a spin system plp^ ] . In this way, a 
correspondence is established between the statistics of branched, polydisperse, macromolecules and a multi-component 
field theory. This field theory reduces to the one-state limit of the Potts model under circumstances appropriate for 
the transition to a physical gel (i.e. a state in which one is certain to find a reversibly-bonded, infinite, branched 
macromolecule) p^ . 

An essential ingredient of the approaches discussed in the previous paragraph is the Potts model in its one-state 
limit — a representation of percolation . It is therefore worth considering similarities and d ifferences between 



the minimal field theory of the vulcanization transition focused on in the present Paper, Eq. (2.7b), and the minimal 
field theory for the Potts model. The minimal field theory for the Potts model is the n ^ limit of the cubic 
(n -|- l)-component field theory, the Landau- Wilson Hamiltonian for which is 

r " / 1 1 \ " 

/ rf'^(E U^^' + s'^^-n-^^'^ E ^S.V'o^V'.), (6.2) 



where r controls the bond-occupation probability (and hence the percolation transition), ' is the nonlinear coupling 
strength, and A^,^^ is the "Potts tensor" (which controls the internal symmetry of the theory; for a discussion of this 
theory see, e.g.. Sec. 2.7 of Ref. 

How does this Potts field theory compare the vulcanization field theory that we have been analyzing in the present 
Paper? The Potts field theory has, a cubic interaction, as does the vulcanization field theory, and therefore its upper 
critical dimension is also six. If we examine the RG analysis of the Potts field theory (in an expansion around six 
dimensions) p6[ | we see that, at the one-loop level, diagrams identical in form (i.e. those shown in Fig. |^) enter the 
renormalization of the various vertices. Moreover, in the n —^ limit the RG flow equations for the two theories turn 
out to be identical. This striking result is connected to the following observations: 

(i) In Potts case, aside from the d-dimensional integrals corresponding to the diagrams, the coefficients in the flow 
equations are determined by the contractions of Potts-tensor indices associated with each cubic vertex, these contrac- 
tions being the origin of the n-dependence of the coefficients in the flow equations. 

(ii) In the vulcanization case, the diagrams intrinsically correspond to (n -I- l)ci-dimensional integrals but, due 
to the constraints on the summations over wavevectors, these diagrams produce (n -I- l)(i-dimensional integrals 
(which smoothly reduces to d-dimensional integrals in the n limit), together with d-dimensional integrals [see 
Eqs. (^,|C|)]. 



(iii) Despite the explicit differences in the forms of the two theories, it turns out that, in the n limit, the integrals 
and the combinatorics conspire to produce precisely the same flow equations. In some delicate way, which we do not 
fully understand, the constraints on the wavevector summations in the vulcanization theory play a similar role to the 
field-index contractions in the Potts theory. 

Having discussed the similarities of the Potts and vulcanization approaches, let us now catalogue the many distinc- 
tions between them: 

(i) The Potts field theory has a multiplet of n real fields on d-dimensional space; the vulcanization field theory has a 
real singlet field living on {n + l)-fold replicated d-dimensional space. 

(ii) The Potts field theory represents a setting involving a single ensemble p^ , the ensemble of percolation configu- 
rations, whereas the vulcanization field theory describes a physical problem in which two distinct ensembles (thermal 
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and disorder) play essential roles. As such, the vulcanization field theory is capable of providing a unified theory not 
only of the transition but also of the structure, correlations and (e.g. elastic) response of the emerging amorphous 
solid state. This is already manifested at the mean-field level, inasmuch as the vulcanization field theory presents an 
order parameter that is far richer in its physical content that the one presented by the Potts model. 

(iii) The entire symmetry structures possessed by the percolation and vulcanization field theories are quite different. 
The Potts field theory has translational and rotational invariancc (in unrcplicated space), along with the discrete 
symmetry of (n + l)-fold permutations of the field components. The vulcanization field theory has the symmetries 
of the independent translations and rotations of the (n + 1) replicas of space, along with the discrete symmetry of 
{n + l)-fold permutations amongst the replicas. 

(iv) The nature of the spontaneous symmetry breaking at the percolation and vulcanization phase transitions is dis- 
tinct. The percolation transition (in its Potts representation) involves the spontaneous breaking of the (n — > limit) 
of a discrete {n + l)-fold permutation symmetry. By contrast, the vulcanization transition involves the spontaneous 
breaking of the (n — > limit of the) continuous symmetry of relative translations and rotations of the n + 1 replicas; 
the permutation symmetry remains intact in the amorphous solid state, as does the symmetry of common translations 
and rotations of replicated space. Thus, the vulcanization transition is associated with the appearance of low-energy, 
long- wavelength, Goldstone-type excitations |2^, which we expect to lead to the restoration of the broken continuous 
symmetry in and below a lower critical dimension of two. By contrast, fiuctuations destroy the percolation transition 
only at and below the the lower critical dimension of unity. 

Whilst there are these apparent distinctions between the percolation and vulcanization approaches, especially in 
low dimensions, there is also evidence in favor of some sort of sharp correspondence between the physics of percolation 
and vulcanization coming from the computation of critical exponents near the upper critical dimension. This apparent 
dichotomy can, however, be reconciled if we carefully delineate between three logically distinct physical properties 
pertaining to RCMSs and other randomly constrained systems: 

(i) macroscopic network formation (by which we mean that constraints are present in sufficient density to connect a 
nonzero fraction of the constituents into a giant random molecule) ; 

(ii) random localization (by which we mean the change in thermal motion of a nonzero fraction of the constituents 
from wandering throughout the container to fluctuating only over finite distances from their random mean positions); 
and 

(iii) the acquisition of rigidity (by which we mean the emergence of a nonzero static shear modulus) . 

Within mean-field theory (and hence above six spatial dimensions), these three properties go hand in hand, emerging 
simultaneously at the phase transition. At and below six dimensions they appear to continue to go hand in hand 
(although, strictly speaking, we have not yet investigated the issue of the acquisition of rigidity beyond mean-field 
theory) until one reaches two dimensions where we believe this broad picture will change (as we shall discuss shortly) . 
Thus, it appears that, within the limited sphere of issues concerning amorphous solidification that percolation-based 
approaches are capable of addressing, such approaches do not lead one astray. In other words, the superposition of 
thermal fluctuations on the positions of the constituents of the macroscopic network that emerges as the constraint 
density is increased towards the phase transition does not lead to any changes in the critical exponents governing 
percolation-type quantities: disorder fluctuations appear to play a more important role than do thermal fluctuations, 
as far as the percolative aspects of the critical phenomenon are concerned. 

This brings up the interesting issue of the nature of the vulcanization transition and its relationship with the 
percolation transition as the dimensionality of space is reduced to the neighborhood of two spatial dimensions, two 
being the lower critical dimension of the vulcanization transition. (The ideas reported in this paragraph result from 
an ongoing collaboration with H. E. Castillo [Q.) Indeed, the case of two dimensions is especially fascinating in view 
of the fact that there is a conventional percolation transition in two dimensions, whereas the thermal fluctuations 
are expected to be sufficiently prominent to destablize the amorphous solid phase, in which case the macroscopic 
network formation no longer occurs simultaneously with the random localization of constituents of the network. It is 
tempting to speculate that in two dimensions an anomalous type of vulcanization transition (not accompanied by 
true localization) continues to happen simultaneously with percolation transition. As the constraint density is tuned 
from below to above criticality, the amorphous solidification order parameter would remain zero, whereas the order- 
parameter correlations would change from decaying exponentially to decaying algebraically with distance. One might 
say that (constraint-density controlled) cluster fragmentation (rather than the thermal excitation of lattice defects, 
as in regular two-dimensional melting) would be mediating the melting transition. If this scenario should happen to 
be borne out then, at sufficiently high crosslink densities one would have a quasi-amorphous solid state — the random 
analogue of a two-dimensional solid |48| — exhibiting quasi-long-range positional order but of a random rather than 
regular type. By implementing these ideas via an effective field theory that describes low-energy excitations of the 
amorphous solid state, we hope to construct a picture of the vulcanization transition and the emergent rigid state in 
and near two spatial dimensions. Such a development would complement the approach to the vulcanization transition, 
presented here, which is based on expanding about the upper critical dimension, e.g., by providing access to critical 
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exponents via an expansion about two rather than six dimensions. 
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APPENDIX A: INVERSE SUSCEPTIBILITY AND GINZBURG CRITERION 

In order to calculate the one-loop correction to the 2-point vertex function r(^'(fc), we first calculate the self- 
energy (i.e. the sum of all two-point one-particle- irreducible amputated diagrams), in terms of which r(^-'(A;) = 
Go(fc)~^ — S„(fc)|^^p. To one-loop order, Sn(fc) is given by the amputated diagram shown in Fig. ^, 



S„(fc) = J2 Go(fci) Go(fci - k). 



(Al) 



fcl eHRS 
(k-hi GHRS) 



Let us emphasize the meaning of the notation: one is directed to sum over all replicated wavevectors fci G HRS subject 
to the constraint that fc — fci e HRS; one should also bear in mind the fact that the external wavevector fc lies in the 
HRS. This constrained summation can be expressed in terms of several unconstrained summations (for cases in which 
k has nonzero entries in at least three replicas, i.e., lies in the 3+RS) as 



J2 o(fci) = Eo(fci)-EEo(fci) 



(fc-felEHRS) 



ki — pe" 



nO(fci) 



fei=0 



a— p 



ki — pe'^ -\-k 



nO{ki) 



ki — k 



(A2) 



for any O(fci). Here, {e"}2=o is the collection of unit vectors in replicated space, so that, e.g., a generic vector p can 
be expressed as J2a=o P"^". When k belongs to the 2RS [e.g. k — (1^, P, 0, • • • , 0)] there is a slight modification of 
Eq. ([A^) and, instead, we have 



E o(fci) = Eo(fci)-EEo(^i) 

fei "=o P 



ki eHRS 
(k-k-i EHRS) 



nO(fci) 



ki—pe° 



ki —k 



nO{ki 



fci=0 



EE 0(^0 



O(fci) 



a— p 

O(fci) 



ki — pe" +/c 



(A3) 



For the moment, let us focus on the case of fc € 3^RS. By making use of Eq. (A2), and subsequently transforming 
each unconstrained summation into an integral, we obtain 

E„(fc) = 18g2 (u"+i j d(«+i)'ifci Go(fci) Go(fci ~k)^2j2vj d^pGoipe") Go(pe" - k) + 2n Go(6) Go(fc)) . (A4) 



The limit of the validity of the Landau theory (i.e. the tree-level approximation) can be ascertained by enquiring 
when the loop corrections to the inverse susceptibility become comparable its tree-level value. Thus we take the 



long-wavelength limit of the correction (Al) via a sequence of wavevectors k lying in the HRS, obtaining 
S„(A:)|^^g = 18.g'(u"+i / S^+'^^^k^ Go{kf - 2{n + 1)V f d^pGo{p f + 2nGo{0f^ 



(A5) 



At this stage, the n ^ limit may be taken [the reason for this is discussed in Sec. II D| , shortly after Eq. ( 2.7b )]. 
In addition, the integral over the {n + l)-fold replicated space goes smoothly into an integral over the ordinary 
(i.e. unreplicated) space. Thus, we arrive at 
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S(fc)|^^Q = lim E„(fc)|^^Q - l^g^{y j d^pGoipf - 2V J d'^pCip)^). (A6) 

From this expression, we see an example of what turns out to be a typical effect of the exclusion of the IRS, viz., 
that it reverses the sign relative to the unconstrained version. By collecting this loop correction together with the 
tree-level inverse susceptibility, we arrive at the result that we shall use to establish the Ginzburg criterion: 



We mention, in passing, that when k lies in the 2RS, we need to use Eq. (A3) instead of Eq. ( [A^ ) in evaluating the 
constrained summation. The resulting two extra terms in turn out to be nonextensive and non-divergent at the 
transition, and thus do not change the result for the Ginzburg criterion. (The appearance of non-extensive terms may 
seem strange, but also occurs in the semi-microscopic theory of RCMSs, where the free energy for the saddle point 
value of the order parameter has a non-extensive part; for a discussion of this issue see Sec. 2.6 of Ref. B.) 



APPENDIX B: SUBLEADING ELEMENTS: ADDITIONAL SEMI-MICROSCOPICALLY GENERATED 

FIELDS AND VERTICES 



The inspiration for the minimal model, Eq. (2.7b), discussed in Sec. II D, comes from experience with detailed 
statistical- mechanical investigations of various semi-microscopic models of RCMSs and related systems ||,^,|ll|,|l2| . 
The field theories obtained in these investigations contain additional fields and vertices beyond those featuring in the 
minimal model. Amongst them are: the IRS field [variously denoted as r2(ke") or ri"(k)], which describes density 
fluctuations; various vertices that couple the IRS field to itself and to the HRS field; and quartic or higher-order 
HRS vertices. In the present section we discuss the role of these additional fields and vertices. We shall confine our 
attention to effects that show up at the one-loop level. To avoid confusion we shall, in this section, denote the bare 
HRS and IRS correlators respectively by G™^ and Gj^^. 




FIG. 4. Example of a one-loop correction to the self-energy due to a vertex omitted from the minimal model. 



1. Subleading influences on the higher replica sector 

We begin by considering the possible corrections to the HRS self-energy E„(fc)|^_^Q arising from the additional fields 
and vertices. At the one-loop level, the only contributions arising from an omitted vertex are those associated with 
the quartic vertex, for which there are two situations to consider, depending on whether the loop wavevector lies in 
the IRS or the HRS. Figure || shows the relevant diagram. 

Let us first look at the contribution of this diagram when the loop wavevector lies in the HRS. In this case, evaluating 
the diagram involves the constrained summation: 

n 

E Gr'Ck) = E GT'Ck) -EE GTH^n + "G™^(6) 

feeHRS k k 

= E GT'Ck) - (r. + 1) E G^'W + nCrHO), (Bl) 

k k 
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which vanishes in the n — > hmit. 

Let us now look at the contribution of this diagram when the loop wavevector hes in the IRS. In this case, no 
critical bare correlators feature, so that the resulting contribution to I]„(A;)|j,^g is finite. There are also contributions 
to S„(A;)| arising from one- loop diagrams involving two cubic vertices, in which either one or b oth loop-wavevectors 



lie in the IRS. None of these contributions alters the Ginzburg criterion established in Sec. V A 



2. Absence of feedback of critical fluctuations on the density-density correlator 



As we have discussed in Sec. HI, the IRS field il(ke"), which describes density fluctuations, remains "massive" at 
the vulcanization transition (i.e. the coefficient of the term quadratic in this field remains positive at the transition), 
and the corresponding bare correlator is nonsingular at the vulcanization transition. We now examine the effects 
of HRS critical fluctuations on the correlator of the IRS field. We approach this issue by studying those one-loop 
diagrams for the IRS self-energy in which at least one internal wavevector lies in the HRS; there are three types of 
contribution to consider: 

(i) There is the contribution associated with the diagram shown in Fig. ^ but with the external wavevectors now lying 



in the IRS. By the same reasoning that we applied to Eq. (Bl), this contribution vanishes in the n ^ limit, 
(ii) There are the two contributions associated with the type of diagram shown in Fig. |l]. When one of the internal 
wavevectors lies in the IRS and the other lies in the HRS, the contribution involves a constrained summation over 
k with k € HRS and (pe" — k) G IRS (where k is the loop wavevector and pe" is the external wavevector). In this 
case, the constraints on the summation require that k £ 2RS and k — pe" + le^ , where f3 ^ a and 1 7^ 0. Then, the 
contribution to the IRS self-energy reads 

Y: G^Hk) (pe- - fc) = E E GT'ipe'^ + le^) (-le'') 

= nEGr'(pe"+le'')Gji^S(_i-/3)|^^^ ^ (g2) 
15^0 

which evidently vanishes in the n — s- limit. On the other hand, when both internal wavevectors lie in the HRS, 
the contribution involves the constrained summation over k with k g HRS and (pe" — fc) e HRS. In this case, the 
contribution to the IRS self-energy reads 

E G™S (fc) G™S (pga -k)= (fc) Gr' (Pe" - fc) - E (^) (P^" - fc) oc n , (B3) 

fceHRS fceHRS feeHRS 

(peO-fceHRS) (pc-O-fceiRS) 



which also evidently vanishes in the n ^ limit. [In the last step we have used Eq. (B2), as well the strategy for 



handling constrained summations employed in Eq. (Bl) 



We conclude that, to one-loop order, the IRS self-energy does not acquire any singular contributions due to critical 
fluctuations in the HRS. In this sense, the two sectors are well separated in the neighborhood of the vulcanization 
transition. However, it is straightforward to show [^9| that there are IRS correlators, such as those involving four 
IRS fields but only two replica indices, which do become long ranged at the vulcanization transition and are thus 
capable of signaling the transition. 

APPENDIX C: DERIVATION OF FLOW EQUATIONS WITHIN THE EPSILON EXPANSION 

1. Implementation of the momentum-shell RG 

The first step in the momentum-shell RG approach that we are adopting is to integrate out Fourier components of 
the field ^l{k) having wavevectors in the shell A/b < \k\ < A. To do this, we define il< and fl^ , respectively the long 
and short wavelength components of Q{k), by 

17<(fc) = (0'. forA/6<|fc|<A; ^ 
^ ^ \n{k), for < |fc| < A/6; ^ ^ 

^>rk)^lm, forA/6<|fc|<A; 

^ ' lO, for < |fc| < A/b. ^ ' 
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Then, by exchanging il{k) for ^-^{k) and r2<(fc) in Eq. (5.6) we can re-express the effective Hamiltonian as 



Sn{m) = s4{n<}) + N 

fceHRS 

V{m) ^Ng ^ 



+ -m\n>(k)\'~v{{n}), 



(C2a) 



fci,fe2,fc3eHRS 

+ 3n<{ki) n>{k2) ^^ih) + 3r2<(fci) n<{k2) i^^ih) 



(C2b) 



Now, focusing on the partition function, we integrate out the aforementioned short-wavelength field components fi^ 
in the context of a cumulant expansion in V. Thus, Eq. (2.7a) becomes 



[Z" 



^ V'n< exp ( - 5<'^*^), 
5<'<=«({17<}) ^ 5„({17<}) - ln(cxp V^)^ , 

V^'n^ exp (-7V^. 



(C3a) 
(C3b) 



In ^ exp y ) = In 



fceHRS 



-\k\')\n>{k)\ expF 



J V^n> exp(-7V^. 



feeHRS 



-\k\^)\n>{k)\ 



- (^)> + ^((^')> - (^)>) + i((^')> - 3(^)>(^^)> + 2(y)t) + OiV^ (C3c) 

Note that we have not explicitly given the factor associated with Gaussian fluctuations in the wavevector shell 
because it is nonsingular and, therefore, does not contribute to the quantities that we are focusing on, viz., the RG 
flow equations. 

Next, we calculate 5^'°* to the one-loop level by computing the cumulant expansion to 0{V^) and discarding 
operators that are irrelevant in the vicinity of d = 6. This amounts to retaining only terms of the form of those 
present in the original minimal model, and thus we are in a position to begin the task of recasting the resulting 
theory in its origi nal form. The terms that must be considered correspond to the diagrams shown in Fig. and are 
computed in Sec. C2, When included, they produce the following intermediate form for the effective coarse-grained 
Hamiltonian: 



+k2+k3,0 ■ 



(C4) 



feGHRS 



/ci,/c2,fc3GHRS 



where the functions /2 and /s can be found in Sec. C2. In fact, only their long wavelength parts are needed, i.e., we 
shall only need the constants /2°\ /2^'' and f^^^ in the Taylor expansions 



/2(fc)«/f + i/2^'^|fcp + 0(fc^), 



(C5a) 
(C5b) 



/3(fci, k, k) ~ /i"^ + oCklkl kl k ■ k, k ■ k, k ■ k)- 

The next step is to rescale il^ and k via 

D.<{k ^ zn'{k'), (C6a) 

= (C6b) 

The reca sting of the theory in its original form also involves the restoration of the wavevector lattice, as discussed in 
Sec. VB. Having made this restoration, we arrive at recursion relations for r and g, along with the condition that 
the coefficient of the gradient term be restored to its original value: 

r' = (r + fi''^/N)zH-^^+^'>-', (C7a) 
g' = (.9 + /f /iV)z36-2(n+i)<i, (C7b) 

1 = (1 - f!^'^/N)z%-'^''+^^''-\ (C7c) 

The computation of the coefficients in the recursion relations simplifies under the convenient choice of 6 = 1 + x 
with X positive and very small, because it allows the approximation of the shell integrals by the product of end-point 
values of the integrands and the shell volumes. Thu s, we arrive at the differential RG rec ursion relations (i.e. flow 
equations) given in the main text in Eqs. (5.8a) and (5.8b), along with the coefficients (^). 
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2. Evaluation of two diagrams 



The renormalizations of t and the gradient term acquire a nontrivial contribution associated with diagram (a) of 
Fig. ^, which determines f2{k) in Eq. (C4). Thus, including the symmetry factor of the diagram, we need to evaluate 

/2(fc) = V E Go(fci)Go(fci-fc). (C8) 

ki GHRS 

We have encountered this kind of constrained summation in App. and we use the recipe given there, together with 
the facts that the external wavevector satisfies |fc| G (0, A/6) whereas the internal wavevectors satisfy \ki \ £ (A/6, A) 
and |fci — fc| G (A/6, A). In practice, we are concerned with the small-fc behavior of f2{k), in which case the latter two 
constraints are equivalent (the difference in their effects being sub-dominant). Thus, by invoking Eq. (A2) we arrive 
at 



/ 



yn+l 



{2n) 



(n+l)d 



d(^+i)dk^ Go(fci)Go(^i -fc) 



-2E 



A/b<|fc|<A 
V 



(2^) 



d''pG'o(pe")Go(pe"-A:) 



(C9) 



A/6<|p|<A 

Then, by expanding for small k and taking the n limit, we obtain 

k'^-^dk 



AO) _ " „2 
J2 — 



4^ {2irY (-r + k^/2Y 



dk 



2 k'^+^dk 



C^ttY V JA/b i-r + ky2r d J^/, (-r + P/2)4 



(ClOa) 
(ClOb) 



where Sd is the surface area of a d-dimensional sphere of unit radius. 

The renormalization of g acquires a nontrivial contribution associated with diagram (b) of Fig. |^, which determines 
f^{ki, ^2, fcs) in Eq. (p4|). Thus, including the symmetry factor of the diagram, we need to evaluate 



hikMM) ^ j^iigf J2 GoCk)GoCk + k2)GoCk-k) 



(Cll) 



fcGHRS 
(fe + fc2eHRS) 
(fc-fcieHRS) 



This constrained sum is similar to the one analyzed in the context of Eq. (A2), but is more lengthy, yielding 

yn+l 



/3(fci,fc2,fc3) =36g^ 



P^^in+l)d 
n 

-E 



A/b<|fe|<A 
V 



a=0 ^ ' 



^(n+l)d^ Go(fc) Go(fc + fc2) Go(fc - fcl) 

d'^pGo(fc) Go(fc + fc2) Go(A: - 



A/f)<|p|<A 



E 



V 



(27r)'* 



d^'pGo(fc) Go(fc + fc2) Go(A: - fci)|^.^ 



A/b<|p|<A 



V 



' ^ {2t:Y 

a=0 ^ ' 



d'^pGo(fc) Go(fc + fc2) G^(k - fci)|fc^p,„+fe^ 



(C12) 



A/b<|p|<A 

In fact, what we need is the n — > limit of /3(0, 0, 0), which is given by 
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/f = -9g3^^ r . + (C13) 



(27r)''7A/6 (-r + fcV2) 



It is worth emphasizing that, once again, the essential consequence of the exclusion of the IRS from the theory. 
Without it, even signs of all three co^ 
theory would be completely changed. 



Without it, even signs of all three coefRcients, f2°\ /j^^ and /a, would be reversed, and the fixed-point structure of 
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